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WIENER HOPF FACTORIZATION AND DISTRIBUTION OF 
EXTREMA FOR A FAMILY OF LEVY PROCESSES 

By Alexey Kuznetsov 1 

York University 

In this paper we introduce a ten-parameter family of Levy pro- 
cesses for which we obtain Wiener-Hopf factors and distribution of 
the supremum process in semi-explicit form. This family allows an 
arbitrary behavior of small jumps and includes processes similar to 
the generalized tempered stable, KoBoL and CGMY processes. An- 
alytically it is characterized by the property that the characteristic 
exponent is a meromorphic function, expressed in terms of beta and 
digamma functions. We prove that the Wiener-Hopf factors can be 
expressed as infinite products over roots of a certain transcendental 
equation, and the density of the supremum process can be computed 
as an exponentially converging infinite series. In several special cases 
when the roots can be found analytically, we are able to identify 
the Wiener-Hopf factors and distribution of the supremum in closed 
form. In the general case we prove that all the roots are real and 
simple, and we provide localization results and asymptotic formu- 
las which allow an efficient numerical evaluation. We also derive a 
convergence acceleration algorithm for infinite products and a sim- 
ple and efficient procedure to compute the Wiener-Hopf factors for 
complex values of parameters. As a numerical example we discuss 
computation of the density of the supremum process. 

1. Introduction. Wiener-Hopf factorization is a powerful tool in the 
study of various functionals of a Levy process, such as extrema of the process, 
first passage time and the overshoot, the last time the extrema was achieved, 
etc. These results are very important from the theoretical point of view; for 
example, they can be used to prove general theorems about short/long time 
behavior (see [7, 18, 21] and [29]). However, in recent years, there has also 
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been a growing interest in applications of Wiener-Hopf factorization, for 
example, in Insurance Mathematics and the classical ruin problem (see [2]) 
and in Mathematical Finance, where the above-mentioned functionals are 
being used to describe the payoff of a contract and the corresponding prob- 
ability distribution is used to compute its price (see [3, 10, 26] and [30] and 
the references therein). 

Let us summarize one of the most important results from Wiener-Hopf 
factorization. Assume that Xt is a one-dimensional real- valued Levy process 
started from Xq = and defined by a triple (p,, a, u), where [i E K specifies 
the linear component, a > is the volatility of the Gaussian component and 
v(dx) is the Levy measure satisfying f^mm{l,x 2 }i/(dx) < oo. The charac- 
teristic exponent ^(z) is defined by 

E[e izXt ] = e~ t ^ z \ zeR, 

and the Levy-Khintchine representation (see [7]) tells us that ^(z) can be 
expressed in terms of the generating triple (/x, a, v) as follows: 

(1) V(z) = \a 2 z 2 -ifiz- [ {e izx - 1 - izh{x))u{dx). 

Here h(x) is the cut-off function, which in general can be taken to be equal to 
ccl{| a .| < ij.; however, in this paper we will use h(x) = (Section 2) or h(x) = x 
(Sections 3 and 4). 

We define extrema processes 

S t = sup{X s : < s < t}, I t = inf{X s : < s < t} 

introduce an exponential random variable r = r(q) with parameter q > 0, 
which is independent of the process Xt, and use the following notation for 
characteristic functions of S T and I T : 

<f>+(z) = E[e izS -M], <j)-(z)=E[e izI ^]. 

The Wiener-Hopf factorization states that the random variables S T and 
X T — S T are independent, random variables I T and X r — S T have the same 
distribution; thus for zGKwe have 



izX T ] 



q + ^(z) 
(2) 

= E[e i2 ^]E[e i2(XT - 5 ^] = <f>+(z)<l>-(z). 

Moreover, random variable S T (J T ) is infinitely divisible, positive (negative) 
and has no linear component in the Levy-Khintchine representation (1). 
There also exist several integral representations for in terms ofF(Xt E dx) 
(see [7, 18, 21] and [29]) or in terms of 9{z) [6, 25]. 
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The integral expressions for the Wiener-Hopf factors 4>q are quite compli- 
cated; however, in the case of stable process it is possible to obtain explicit 
formulas for a dense class of parameters (see [17]). It is remarkable that in 
some cases we can compute Wiener-Hopf factors explicitly with the help of 
factorization identity (2). As an example, let us consider the case when the 
Levy measure is of phase- type. Phase- type distribution (see [2]) can be de- 
fined as the distribution of the first passage time of a finite state continuous 
time Markov chain. A Levy process Xt whose jumps are phase-type dis- 
tributed enjoys the following analytical property: its characteristic function 
^(z) is a rational function. Thus function q(q + ^(z)) -1 is also a rational 
function, and therefore it has a finite number of zeros/poles in the complex 
plane C. And here is the main idea: since the random variable S T (I T ) is 
positive (negative) and infinitely divisible, its characteristic function must 
be analytic and have no zeros in C + (C~), where 

C+={zeC:Im(z)>0}, C~ = {z G C : Im(z) < 0}, C ± = C ± UR. 

Thus we can uniquely identify 4>~^(z) [4>q{z)] as a rational function, which 
has value one at z = and whose poles/zeros coincide with poles/zeros of 
qiq + ^iz))" 1 in C" (C+). 

While Levy processes with phase-type jumps are very convenient objects 
to work with and one can implement efficient numerical schemes, there are 
some unresolved difficulties. One of them is that by definition phase-type 
distribution has a smooth density on [0, oo); in particular the density of 
the Levy measure cannot have a singularity at zero. This means that if 
we want to work with a process with infinite activity of jumps, we have 
to approximate its Levy measure by a sequence of phase- type measures, 
but then the degree of rational function ^(z) would go to infinity and the 
above algorithm for computing Wiener-Hopf factors would quickly become 
unfeasible. 

In this paper we address this problem and discuss Wiener-Hopf factor- 
ization for processes whose Levy measure can have a singularity of arbitrary 
order at zero. The main idea is quite simple: if characteristic exponent *&(z) 
is meromorphic in C and if we have sufficient information about zeros/poles 
of q + ^f(z), we can still use factorization identity (2) essentially in the same 
way as in the case of phase- type distributed jumps, except that all the finite 
products will be replaced by infinite products, and we have to be careful 
with the convergence issues. The main analytical tools will be asymptotic 
expansion of solutions to q + ^f(z) = and asymptotic results for infinite 
products. 

The paper is organized as follows: in Section 2 we introduce a simple 
example of a compound Poisson process, whose Levy measure has a den- 
sity given by v(x) = e aa; sech(x). We obtain closed form expressions for the 
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Wiener-Hopf factors and density of S T . Also, in this simple case we intro- 
duce many ideas and tools which will be used in other sections. In Section 3 
we introduce a Levy process Xt with jumps of infinite variation and the den- 
sity of the Levy measure v{x) = e ax sinh(:r/2) -2 . This process is a member 
of the general /3-family defined later in Section 4; however, it is quite unique 
because its characteristic exponent *&(z) is expressed in terms of simpler 
functions, and thus all the formulas are easier and stronger results can be 
proved. In this section we derive the localization results and asymptotic ex- 
pansion for the solutions of q + *$>(iz) = 0, prove that all of them are real and 
simple, obtain explicit formulas for sums of inverse powers of these solutions 
and finally obtain semi-explicit formulas for Wiener-Hopf factors and distri- 
bution of supremum S T . In Section 4 we define the ten-parameter /3-family 
of Levy processes and derive formulas for characteristic exponent and prove 
results similar to the ones in Section 3. Section 5 deals with numerical issues: 
we discuss acceleration of convergence of infinite products and introduce an 
efficient method to compute roots of q + ^f(z) for q complex. As an example 
we compute the distribution of the supremum process St- 

2. A compound Poisson process. In this section we study a compound 
Poisson process Xt, defined by a Levy measure having density 

v(x) 



cosh(x) 

We take the cut-off function h(x) in (1) to be equal to zero, and thus the 
characteristic exponent of Xt is given by 

TT TT 



(3) Wz) = - / (e lxz - lWx) dx - 

w v ; V ' y 1 cos(vr/2a) cosh(vr/2(z - ia)) 

and the above integral can be computed with the help of formula 3.981.3 in 
[20]. Our main result in this section is the following theorem, which provides 
closed-form expressions for the Wiener-Hopf factors and the distribution of 

Sr. 

Theorem 1. Assume that q>0. Define 

2 ( TT 

rj = — arccos 



tt \q + TTsec(Tr/2a 

(4) r(l/4(l - a))r(l/4(3 - a)) 

Po ~ 



r(l/4fa-a))r(l/4(4-7j-a))- 
Then for Im(z) > (a — n) we have 

r(l/4(i7 -a- iz))r(l/4(4 -rj-a- iz)) 



(5) #(*)=l>o- 



r(l/4(l -a- zz))r(l/4(3 -a- iz)) 
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We have P(5 T = 0) = po, and the density of S T is given by 

^nsr<x) 

dx 

7T C ° V 2 



(6) 



r(l/4(l + ry))r(l/4(3 + ry)) 
r(l/2r ? ) 

r(l/4(5- 77)^(1/4(7-7?)) 



r(l/2(4-77)) 



5 — 77 7 — T] 4 



V -4* 



where 2-^1 (a, b; c; z) is the Gauss hypergeometric function. Ifq = and a < 0, 
equation (4) implies 77 = |a|, and formulas (5) and (6) are still valid. In this 
case the random variable >S r (o) should be interpreted as Soo = sup{X s : s > 0}. 

First we will state and prove the following lemma, which will be used 
repeatedly in this paper. It is a variant of the Wiener-Hopf argument, which 
we have borrowed from the proof of Lemma 45.6 in [29]. 

Lemma 2. Assume we have two functions f + (z) and f~{z), such that 
/^(O) = 1, / (z) are analytic in C^, continuous and have no roots in C 1 * 1 



and z 1 ln(/ (z)) — > as z 
(7) 



q + ^{z) 



00, zeC 1 // 
= f + (z)f-(z), 



z G 



then f ± (z)=^(z). 



PROOF. We define function F(z) as 



F(z) 



f -(*)'■ 



if z e 



if zeC+. 



W3W 

Function F(z) is well defined for z real due to (7) and (2). Using properties of 
<fy^ and f^ we conclude that F(z) is analytic in C + and C - and continuous 
in C, and therefore by analytic continuation (see Theorem 16.8 on page 
323 in [28]) it must be analytic in the entire complex plane. Moreover, by 
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construction function F(z) has no zeros in C, thus its logarithm is also an 
entire function. All that is left to do is to prove that function ]n(F(z)) is 
constant. 

Using integration by parts and formula (1) one could prove the following 
result: if £ is an infinitely divisible positive random variable with no drift and 
^^(z) is its characteristic exponent, then z~ l ^^{z) — > as z — > oo, z 6 C + 
(this statement is similar to Proposition 2 in [7]). Thus 

z _1 ln(^(z)) ->0, z^oo, zeC 1 . 

Since functions also satisfy the above conditions, we find that z~ x ln(F(z)) 
as \z\ — > oo in the entire complex plane. Thus we have an analytic func- 
tion ln(F(z)) which grows slower than \z\ as z — > oo, and therefore we can 
conclude that this function must be constant (a rigorous way to prove this 
is to apply Cauchy's estimates, see Proposition 2.14 on page 73 in [16]). The 
value of this constant is easily seen to be zero, since /^(O) = ^(0) = 1. □ 

Proof of Theorem 1. Using expression (3) for ty(z) we find that 
function q(q + ^(z)) _1 has simple zeros at + a + An) , i(3 + a + An)} and 
simple poles at {i(a + n + An), i(a — w + An)}, where n € Z and n is defined 
by (4). Next we check that \a\ < rj < 1 and define function f + (z) as product 
over all zeros/poles lying in C~ 



( 8 ) = n a 



(1 - iz/(An + 1 - a))(l - iz/(An + 3 - a)) 



iz/{An + n — a))(l — iz/(An + 4 — n — a)) 

and similarly f~{z) as product over zeros/poles in C + . It is easy to see 
that the product converges uniformly on compact subsets of C \ iM since 
each term is 1 + (9(n -2 ) (see Corollary 5.6 on page 166 in [16] for sufficient 
conditions for the absolute convergence of infinite products). The fact that 
f + (z) is equal to the right-hand side of formula (5) can be seen by applying 
the following result from [19]: 



m IIi 



x/(n + a) r(a)r(6 + a;) 



, x/(n + b) T(b)T(a + x)- 

The formula for f~(z) is identical to (5) with (z,a) replaced by {—z, —a). 

Now we will prove that f ± {z) = (j)^{z). First, using the reflection formula 
for the gamma function (formula 8.334.3 in [20]), one can check that for zGl 
functions f^(z) satisfy factorization identity (7). Next, using the following 
asymptotic expression (formula 6.1.47 in [1]): 



(10) 
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we conclude that z~ x ln(/ (2)) — > as z — > 00, z G C^, and thus all the 
conditions of Lemma 2 are satisfied, and we conclude that f^{z) = 4>f(z). 

To derive formula (6) for the density of S T we use equations (5) and (10) 
to find that E[e~^ s '' r ] = 4>q{i() — >• Po as ( — > 00, where po is given by (4). 
This implies that distribution of S T has an atom at x = (which should not 
be surprising since X t is a compound Poisson process), and P(S T = 0) =po- 
The density of S T can be computed by the inverse Fourier transform 

d -P(S T < x) 



dx 



1 

Po_ 
2vr 



[^(z)-p ]e-^dz 
r(l/4(7? -a- iz))r(l/4(4 - r] - a - iz)) 



1 



' dz. 



r(l/4(l - a - iz))r(l/4(3 - a - iz)) 

Formula (6) is obtained from the above expression by replacing the contour 
of integration by ic + M, letting c — > —00 and evaluating the residues at 
z G {—i{An + T) — a), — i{An + 4 — rj — a)} for n > 0. Evaluating the residues 
can be made easier by using the reflection formula for the gamma function. 
□ 



Remark 1. There are other examples of Levy measures v{x)dx, which 
have finite total mass (and thus can define a process with a finite intensity of 
jumps), and for which the characteristic exponent is a simple meromorphic 
function. These are two examples based on theta functions (see Section 8.18 
in [20] for definition and properties of theta functions): 



vi(x) 



e- ax e 2 (0,e 



e- ax 3 (O,e 



n>0 



-(n+l/2) 2 x 



n>0 



These two jump densities are defined on x > 0, they decay exponentially as 
x — > +00 and behave as x~ l l 2 as x — > + ; thus the total mass is finite. The 
Fourier transform of these functions can be computed using formulas 6.162 
in [20] 



*'v\(x) dx 



z V2{x) dx 



7T 



\Ja — iz 

7T 



tanh(7T\/a — iz) , 
coth(-7r\/a — iz). 



/o y/oT— iz 

Unfortunately equation q + ^f(z) = cannot be solved explicitly which im- 
plies that we cannot obtain closed form results as in Theorem 1; however, 
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these processes could be treated using methods presented in the next sec- 
tions. 



3. A process with jumps of infinite variation. In this section we study a 
Levy process X t , defined by a triple (fi, a, u), where the density of the Levy 
measure is given by 



[sinh(>/2)] 2 

with |q| < 1 (it is a Levy measure of a Lamperti-stable process with char- 
acteristics (1,1 + a, 1 — a), see [22]). The jump part of X t is similar to 
the normal inverse Gaussian process (see [4, 15]), as it is also a process of 
infinite variation, the jump measure decays exponentially as \x\ — > oo and 
has a 0(x~ 2 ) singularity at x = 0. Note that since the Levy measure has 
exponential tails we can take the cut-off function h(x) = x in (1). 

By definition process Xt has three parameters. However, if we want to 
achieve greater generality for modeling purposes, we could introduce two 
additional scaling parameters a and b > and define a process Y t = aX^ t , 
thus obtaining a five parameter family of Levy processes. 



Proposition 3. The characteristic exponent of Xt is given by 



(11) ^>(z) = \o- 2 z 2 + ipz + 4tt(z - ia) coth(7r(2i - ia)) - 4j 



where 



I \ A 1 47(7 - 1) 

7 = 7racot(7ra), p = Air a + 



a 



Proof. We start with the series representation valid for x > 0, 

-2 



(12) 



sinhf — 
V2 



"(1 -e~ x ) 



x\2 



ne 



n>l 



which can be easily obtained using binomial series or by taking derivative of 
a geometric series. The infinite series in (12) converges uniformly on (e, 00) 
for every e > 0, thus 



1 — izx 



n>l 

4 E 

n>l 



sinh(x/2) 2 
n n 



dx 



111 z 



n — a — iz n — a (n — a) 2 
a + iz a + iz iaz 



n — a — iz n — a (n — a) s 
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The integral in the Levy-Khintchine representation (1) for *&(z) can now be 
computed as 



oo 



(e zzx - 1 - izx) . —, ,„ dx + / (e IZX - 1 + izx) . - . . = dx 
smn(x) z Jo svnhyxy 

= 8(a + izf — 5 — rr - 8(a + izJaV^ — — ^ — , 

^ n 2 - (a + zz) 2 ^n 2 -a 2 
n>i v ; n>l 

-4zaz> -k H . 

^— ' (n — ar a 

To complete the proof we need to use the following well-known series expan- 
sions (see formulas 1.421.4 and 1.422.4 in [20]) 

, . , 1 1 

COth(7TXj = 1 > ;r, 

7TX TT X z + U z 

n>l 

1 + cot(7rx) 2 = cosec(7rx) 2 = — =■ 777. 

7r 2Z ^(n-x) 2 n 

Note that it is impossible to find solutions to q + ^{z) = explicitly in the 
general case, even though the characteristic exponent ^f(z) is quite simple. 
It is remarkable that in some special cases, when a = and parameters /j, 
a and q satisfy certain conditions, we can still obtain closed- form results. 
Below we present just one example of this type. 

Proposition 4. Assume that a = a = 0. Define 

(13) 7? = — arccot ( 

ir \4ir 

Then Wiener-Hopf factor <j>g(z) can be computed in closed form when q = A, 
The density of S r ^ is given by 

A P(Sr(4) < x) = !i^M^_i r ,. 

dx w 7T 



The proof of Proposition 4 is identical to the proof of Theorem 1. 

Remark 2. This result is very similar to Proposition 1 in [11], where the 
authors are able to compute the law of I T (q) m closed form only for a single 
value of q, and this law is essentially identical to the distribution of <S r (4). 
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This coincidence seems to be rather surprising, since these propositions study 
different processes: our Proposition 4 is concerned with a Lamperti-stable 
process having characteristics (1,1,1) (see [22]) and completely arbitrary 
drift, while Proposition 1 in [11] studies a Lamperti-stable process with 
characteristics (a,l,a) but with no freedom in specifying the drift, which 
must be uniquely expressed in terms of parameters of the Levy measure. 

The following theorem is one of the main results in this section. It de- 
scribes various properties of solutions to equation q + ^f(z) = 0, which will 
be used later to compute Wiener-Hopf factors and the distribution of the 
supremum process. 

Theorem 5. Assume that q > and that ^(z) is given by (11). 

(i) Equation q + = has infinitely many solutions, all of which 
are real and simple. They are located as follows: 

Co" € (a -1,0), 

C + G(0,a + 1), 

(15) 

Cn £ in + a, n + a + 1), n>l 
Cn £ in + a — l,n + a), n < — 1 . 

(ii) If a we have as n — > ±oo 



Cn = in + a) H — p(n + a)" 1 



(16) 



a 2 



, ^ + a ) in + ay 2 + 0(n" 3 ). 



(17) 



(hi) If a = we /iaue as n — >• ±oo 

Cn+5 = (n + a + w ) + c (n + a + wo)" 1 

CO 



■iAj-q-A7T z co)in + a + u} ) ^ + 0(n" 



where 



4(4 7 - g + ap) 1 / p \ 

and <5 G { — 1,0,1} depending on the signs of n and p. 

(iv) Function g(a + ^(z)) -1 can &e factorized as follows: 

(18) Q = 1 TT 1 + iz/jn + a) 

{ ' Q + *{') (l + iz/C+Xl + iz/Q)^ l+iz/( n ' 

where the infinite product converges uniformly on the compact subsets of the 
complex plane excluding zeros/poles ofq + ^(z). 
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First we need to prove the following technical result. 



Lemma 6. Assume that a and B are not equal to a negative integer, and 



0(n S1 ) for some e\ > as n — > oo. Then 



(19) 



n r 

n>0 



1 + z/(n + a) 
+ z/(n + B + b n ) 



Cz* 



as z 



oo, |arg(»| <7r-e 2 <vr, where C = pgj ]ln>o( 1 + 



Proof. First we have to justify absolute convergence of infinite prod- 
ucts. A product in the left-hand side of (19) converges since each term is 
1 + 0(n~ 2 ) and the infinite product in the definition of constant C converges 
since each term is 1 + 0(n~ 1_£l ) (see Corollary 5.6 on page 166 in [16] for 
sufficient conditions for the absolute convergence of infinite products). Thus 
we can rewrite the left-hand side of (19) as 



(20) 



III 

n>0 



1 + zj(n + a) 
+ z/{n + B + b n ) 



IIi 



1 + z/(n + a) 



yr l+z/(n + P) 

_ + z/{n + B) LL i + z/in + P + br, 

n>0 n>0 



c 



T(B + z) 



n 



z + n + B 



T(a + z) L *-z + n + B + b n 

n>0 



The ratio of gamma functions gives us the leading asymptotic term z^~ a 
due to (10). Now we need to prove that the last infinite product in (20) 
converges to one as z — > oo, |arg(z)| < ir — e 2 < it. We take the logarithm of 
this product and estimate it as 



n>l 
n>l 



z + n + B 

\b n \ 

\z + n + /3\ 



~ ^ \z + n + B\ 

n>l 1 M ' 
1 



< 



^n £ ^\z + n + B\ 



where we have used the fact that ln(l + x) < x for x > and |6 re | < An £l 
for some A > 0. Since |arg(z)| < it — e 2 < vr we have for z sufficiently large 
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\z + n + (3\ > max{l, |n — \z + Let m = [\z + /3|], where [x] denotes the 
integer part of x. Then 

^ 21 ^ ^ n £ iU + n + B\ < ^ n £ ifm + 1 -n) + ^ 



n>l 



+ n + /3| ^-^ n £l (m + 1 — n) ^— ' n £l (n — m) 

rl n=l y ' n=m+l y ' 



The first series in the right-hand side of (21) converges to zero as m — > 00, 
since 



n £l (m + 1 — n) n £l (m + 1 — n) 

n=l v 7 n=l v ; re=[y^n]+l V 7 



< IvH + m -i/2 y 1 

m + 1 — [\/rn\ ^— ' (m + 1 — n) 

n=[ v / m] + l 



< Lvm = +m" £l / 2 ln(m). 

m + 1 — |y ?nj 

The second series in the right-hand side of (21) can be rewritten as 
Y^=ii n + m)~ £l n~ 1 , and we see that it is a convergent series of positive 
terms, where each term converges to zero as m — > 00. By considering its 
partial sums it is easy to prove that the series itself must converge to zero 
as m — > 00 . □ 



Proof of Theorem 5. The proof consists of three steps. The first step 
is to study solutions to equation q + = 0. We will produce a sequence 

of "obvious" solutions Q n and study their asymptotics as n — > ±00. Note 
that this first step requires quite demanding computations, which can be 
made much more enjoyable if one uses a symbolic computation package. 
The second step is to represent the function q(q + ^(z)) _1 as a general 
infinite product, which includes poles of ^(z), zeros of q + ^(z) (given by 
"obvious" ones {«CoS*Cn} an d possibly some "unaccounted" zeros) and an 
exponential factor. Our main tool will be Hadamard theorem (Theorem 
1, page 26, in [24] or Theorem 3.4, page 289, in [16]). We produce entire 
functions P{z) and Q(z), such that P{z) has zeros at poles of *&{z) and 
q(q + ^(z))~ 1 = P(z)/Q(z). After studying the growth rate oi P(z) and Q(z) 
we apply Hadamard theorem and obtain an infinite product for function 
Q(z) [function P(z) will have an explicit infinite product]. These results give 
us an infinite product for q(q + ^f(z))~ 1 . The last step is to prove the absence 
of exponential factor and "unaccounted" zeros in this infinite product, and 
here the main tool will be asymptotic relation (19) for infinite products 
provided by Lemma 6. 
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First we will prove localization result (15). We use (11) to rewrite equation 
q + y$(i() = as 

(22) 4vr(C - a) cot(7r(C - a)) -(j> + m)C - *7 = ~ K - 1- 

Note that we have separated the jump part of ^(z) on the left-hand side 
and the diffusion part on the right-hand side of (22). See Figure 1, where 
the jump part is represented by black line and diffusion part by grey dotted 
line. 

The left-hand side of (22) is zero at £ = and goes to— oo as ( / a + 1 
or C \ a — 1 (see Figure 1). The right-hand side is negative at £ = and 
continuous everywhere; thus we have at least one solution Q~ G (0, a + 1) and 
at least one solution £ (a — 1,0). In fact it is easy to prove that we have 
exactly one solution on each of these intervals, since 47r(£ — a) cot(7r(£ — a)) 
is a concave function on (a — 1, a + 1), while ^cr 2 £ 2 — ^ — q is convex. 

Next, for n/Owe have 

4"7r(£ — a) cot(7r(C — a)) f~ +oo as £ a — n, £ \ a + n, 
47r(C — a) cot(7r(C — a)) \ — oo as (\a-n,(/a-|-ii, 

thus there must exist at least one zero ( n on each interval (n + a, n + a + 1), 
(n + a — 1, n + a). 



-20 
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Fig. 1. Illustration of the proof of Theorem 5. 
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Next we will prove the asymptotic expansion (16). Since we have assumed 
that a^Owe can rearrange the terms in (22) to obtain 

4(C-a) 



— tan(-7r(C — a)) = — 7 



(23) 



-r 1 



(J 



1-aC 



-l 



1 + 2pa- 2 C- 1 + 0((- 2 ) 



a- 



The main idea in the above calculation is to expand the rational function 
in the Taylor series centered at £ = oo. Now, the right-hand side of (23) is 
small when C, is large, and thus the solution to (23) should be close to the 
solution of tan(7r(£ — a)) = 0, which implies 



(24) 



C = n + a + uj 



and u> = o(l) as n — > oo. Next we expand the right-hand side of (23) in 
powers of w as 

11 , 

— tan(7r(C — a)) = — tan(7ro;) = oj + 0(lu ) 



7T 



7T 



and, using the first two terms of the Maclaurin series for £ 1 in powers of u 

(- 1 = ( n + a + a;) -1 = (n + a)" 1 - w(n + a) -2 + 0(w 2 n -3 ), 

we are able to rewrite (23) as 

8 8 / 2 \ 

u + 0{u 3 ) = — ((n + a)" 1 + w(n + a)~ 2 ) - - + a )(n + a) -2 + 0(n~ 3 ). 
cr z a z \p J 

Asymptotic expansion (16) follows easily from the above formula and (24). 
If a = 0, equation (23) has to be modified as follows: 

- tan(7r(C - a)) = — — = - 

vr pC + ^i-q P 



1 + (4 7 - q)p~ 1 ( 



(25) 



4 4 , 
P P 2 

+ l<!3LZi) (47 _, + £t( , )r2 + 0(r 3 ) , 

where again we have expanded the rational function in the Taylor series 
centered at ( = oo. As before, when £ is large the solution of (22) should be 
close to the solution of 

1 4 

— tan(7r(£ — a)) = -, 

7T p 
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and thus we should expand both sides of (25) in the Taylor series centered 
at the solution to the above equation. We define uj as 

1 /4vr\ 

(26) Q = n + a H — arctan — + uj = n + a + ujq + uj, 

TT V P J 

and again uj = o(l) as n — > oo. To expand the left-hand side of (25) in power 
series in uj we use an addition formula for tan(-) and find that 

1 , „ 1 ( /4vr\ 

L — 



(27) 



— tan(-7r(C — a)) = — tan I arctan) — I + ttuj 

TT TT V V P J 



1 47r/p + tan(7ru;) 1 Att / p + ttuj + 0(uj 3 



{16tt 2 + P 2 ) 



uj H uj 2 



tt 1 — Air/ ptan^TTUj) tt 1 — Att/pttuj + 0(uj 3 ) 

= - + 167T 2 + p 2 )UJ 

p p~ 

+ ^(167r 2 + p 2 V 2 + 0(a; 3 ). 

Again, we use (26) to obtain the Maclaurin series of £ -1 in powers of uj 
C l = (n + a + uj + ujy l 

= (n + a + wo)" 1 - oj(n + a + uj )~ 2 + 0(uj 2 n~ 3 ). 
Using (27) and the above expansion we can rewrite (25) as 

47T 2 
P 

= -4(47 - Q + Oip)((n + a + ujq)' 1 - uj(n + a + ujo)~ 2 ) 

+ 4(47 ~ g) (47 -q + ap)(n + a + w )" 2 + 0(n" 3 ) + 0(w 3 ), 

and from this equation we obtain the second asymptotic expansion (17). 

Now we are ready to prove the factorization identity (18) and the fact 
that all the zeros of q + ^f(iC) are real and simple and that there are no 
other zeros except for the ones described in (15). First we need to find an 
analytic function P(z) such that P (0) = 1 and which has zeros at all poles 
of ^(z) (with the same multiplicity). The choice is rather obvious due to 
(11): 

(28) x smhKz-za)) ^ 

sin(-7ra) z — ia 

By definition, the function 

(29) Q(z)=q~ l (q + *(z))P(z) 



16 



A. KUZNETSOV 



is also analytic in the entire complex plane. 

Next, using the definition of P{z) (28) and Q(z) (29) we check that Q(z) = 
if and only if q + ^f(z) = 0. We have proved already that the zeros of q + 
$?(iC) include CnjC^j however, some of them might have multiplicity greater 
than one, and there also might exist other roots (real and/or complex). Let 
us denote the set of these unaccounted roots (counting with multiplicity) as 
3- Using asymptotic expansions given by equations (23) and (25) one can 
easily prove that 3 is a finite set (possibly empty). 

Using equations (28), (29) and (11) we obtain an explicit formula for Q{z) 
from which it easily follows that Q(z) has order equal to one, which means 
that one is the least lower bound of all 7 > such that Q(z) = 0(exp(|z| 7 )) 
as z — > 00; the rigorous definition can be found in [24], page 4 or Chapter 
11 in [16]. Since Q(z) has order equal to one, we can use the Hadamard 
theorem (see Theorem 1, page 26, in [24] or Theorem 3.4, page 289, in [16]) 
to represent it as an infinite product over its zeros 

Q(z) = exp(ciz) ( 1 + J ( 1 + 



QJ V Co 

xn(^|)n(i + f)e XP (-f 

^e3 V fc/ M>i V W V U 

for some constant c\ G C. As the next step we rearrange the infinite product 
in the above formula and obtain 

I %z \ t %z 
Q(z) = exp(c 2 z) 1^1 + ^) ( 1 + 

(30) 



CoW Co" 



X 



for some other constant c 2 £ C, where the infinite product converges abso- 
lutely since each term is 1 + 0(n~ 2 ) as n — > 00. Using definition of P(z) (28) 
and infinite product representation of trigonometric functions (see formulas 
1.431 in [20]) we find that 



n>l 



IZ \ IZ 



n + a J \ —n + a 



Combining equations (28), (29), (30) and (31) we finally conclude that for 
all z G C 

(32) - = exp(c 2 z) -pr 1 1 r l + iz/(n + a) 

q + *(z) (l + iz/<£)(l+iz/<£)}^l + iz/z k U 1 + WCn 
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First let us prove that c 2 = 0. Denote the left-hand side of (32) as F\(z) 
and right-hand side as F2(z). Since ^f(z) is a characteristic exponent, it must 
be 0(z 2 ) as z — > oo, z E R, thus clearly z~ l Xn(F\(z)) — > as z — > oo, z G R. 
Using Lemma 6 we find that z~ l \n(F 2 {z)) — > C2 as z — > oo z£l, which 
implies that C2 = 0. 

All that is left to do it to prove that 3 is an empty set. The main tool 
is again Lemma 6. Assuming that (7^0 and using asymptotic expansion 
(16) and Lemma 6 we find that the infinite product in (32) converges to a 
constant as z — > oo, z 6 R. Thus function F 2 (z) ~ A 2 z~ 2 ~ M where M is equal 
to the number of elements in the set 3- However, function F\{z) ~ A\z~ 2 as 
z—too, z£R, thus M = and the set 3 must be empty. In the case a ^ 
the proof is identical, except that both F\(z) and F 2 {z) behave like Az~ l , 
which can be established by the asymptotic expression for £ n given in (17) 
and Lemma 6. □ 

Theorem 5 provides us with all the information about the zeros of q + ^{z) 
that we will need later to prove results about Wiener-Hopf factors and 
perform numerical computations. However, we can also compute explicitly 
the sums of inverse powers of zeros. These results can be useful for checking 
the accuracy, but more importantly, for approximating the smallest solutions 

. We assume that a^O and define for m > 



(33) n m = a— 1 + (Co^r™- 1 + (ar^ 1 + E[c m_1 + or 1 ]- 



Asymptotic expansions (16) and (17) guarantee that the series converges 
absolutely for m > 0, thus the sequence {Q m }m>o is correctly defined. 

Lemma 7. The sequence {£l m }m>o can be computed using the following 
recurrence relation: 



n>l 




'm—n 



m>0 



where coefficients {6 n } n >o o re defined as 



b 2n 



(-l)"- 1 ^ 2 "- 1 
(2n)j 



[n(2n - l)aa 2 + ir 2 a(q + 8n) - 2njp] 



(-l) n 7T 2n 

(2n + l)! 



.2 



(34) b 2n+1 



n(2n + 1) 



7r(47r 2 a 2 + 47 2 - ~/q) 



7T 



+ 7r(2n + l)(47 + ap) . 
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Proof. This statement is just an application of the following general 
result. Assume that we have an entire function H(z) which can be expressed 
as an infinite product over the set of its zeros 3 



Taking derivative of \n(H(z)) we find 



-TO— 1 
Z k 



H>{z) = -H{z)Y J {zk-z)- 1 = -Kz)Y J E 

z fe e3 m>o L z fc e3 

and the recurrence relation for J^ 2fc e3 z k~ m ~ 1 * s obtained by expanding H(z) 
and H'{z) as a Maclaurin series, multiplying two series in the right-hand 
side and comparing the coefficients in front of z m . The statement of Lemma 
7 follows by considering an entire function 

H(z) = qir(z — a)Q(iz), 

where Q{z) is defined by (29). Function H(z) has zeros at {a, and 
one can check that the Maclaurin expansion is given by H(z) = Y2n>o bnZ n 
where coefficients b n are defined in (34). □ 

Finally we can state and prove our main results: expressions for Wiener- 
Hopf factors and density of S T . 

Theorem 8. For q>0 

_ 1 -r-r 1 + iz/jn + a) 

(35) 

Infinite products converge uniformly on compact subsets ofC\iW. The den- 
sity of S T is given by 

(36) -f P(S r < x) = -c^Co^ x ~ £ cfaete, 



where 



(37) 



dx 

k<-l 



-q 1 ~ Co /(n + a) 

n<-l 1_ C /Cn 

l-Ck/ik + a) -q l-Ck/(n + a) 
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Proof. Expressions (35) for Wiener-Hopf factors are obtained using 
factorization identity (18) and Lemmas 2 and 6. Expression (36) for the 
density of S T is derived by computing the inverse Fourier transform via 
residues. □ 

Remark 3. Theorem 8 remains true for q = if fi < 0. In this case 
MX% < and S T — > Soo and I T — > — oo as q — > + . From the analytical point of 
view we have CcT < and Co" = [see Figure 1 and (22)]. If /x = 0, then EX\ = 
and the process X t oscillates; thus Soo = 1^ = oo, which is expressed 
analytically by the fact that Co~ = Co~ = 0- 

4. A family of Levy processes. 

Definition 4. We define a /3-family of Levy processes by the generating 
triple ([i,a,v), where the density of the Levy measure is defined as 

(38) v(x) = ci (1 _ e _ hx)Xl I{x>o} + c 2 (1 _ eP2X)M I{x< } 

and parameters satisfy > 0, /3j > 0, c% > and Aj € (0, 3). This Levy mea- 
sure has exponential tails; thus we will use the cut-off function h(x) = x in 
(!)• 

The /3-family is quite rich: in particular, by controlling parameters Aj, 
we can obtain an arbitrary behavior of small jumps, and parameters a,; 
and Pi are responsible for the tails of the Levy measure (which are always 
exponential). Parameters q control the total "intensity" of positive/negative 
jumps. The processes in /3-family are similar to the generalized tempered 
stable processes (see [15]) which were also named KoBoL processes in [9] 
and [10] 

g — OL-\-X gCt—X 

u ( x ) = c + -^r 1 ^} + c - ur^<o} • 

In fact we can obtain the above measure as the limit of Levy measures 
in /3-family. If we set c\ = c + /3 A +, C2 = c_/3 A ~, a\ = a+(3~ l , = q^/?" 1 , 
Pi = P2 = P and let /3 — > + we see that the Levy measure defined in (38) 
will converge to the Levy measure of the generalized tempered stable process. 
Next, when Ai = A2, the processes in /3-family are similar to the tempered 
stable processes (see [5]). If we restrict the parameters even further, c\ = C2, 
Ai = A2 and fi\ = Pi so that the small positive/negative jumps have the same 
behavior, while large jumps which are controlled by may be different, 
and we obtain a process very similar to the CGMY family defined in [13]. 
Finally, if q = 4, /3j = 1/2, \ = 2 and a± = 1 — a and 02 = 1 + a, we obtain 
the process X t discussed in Section 3. 
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If we restrict parameters as a = 0, p\ = /3 2 , Ai = A2 (and n uniquely spec- 
ified in terms of these parameters), then /3-family reduces to a family of 
Lamperti-stable processes, which can be obtained by Lamperti transforma- 
tion from the stable processes conditioned to stay positive (see the original 
paper [23] by Lamperti for the definition of this transformation and its var- 
ious properties). Spectrally one-sided Lamperti-stable processes appeared 
in [8] and [27], and two-sided processes were studied in [11, 12, 14] and 
[22]. Lamperti-stable processes are a particularly interesting subclass of the 
/3-family since they offer many examples of fluctuation identities related 
to Wiener-Hopf factorization which can be computed in closed form (see 
[11, 14] and [22]). 

In the following proposition we derive a formula for the characteristic 
exponent ^f(z) for processes in the /3-family. As we will see, the characteristic 
exponent can be expressed in terms of beta and digamma functions (see 
Chapter 6 in [1] or Section 8.3 in [20]) 

which justifies the name of the family. 



Proposition 9. If \ e (0,3) \ {1,2}, then 



(40) 



where 



a 2 z 2 



ci 



Pi 



IZ 



+ ipz — -r-B cti — — ; 1 — A 



Pi 



§-B{ a 2 + ^;l-A 2 ] +7, 



P 



7 : 
P 



^B(a 1 ;l-A 1 ) + -^B(a 2 ;l-A 2 ) 
Pi P 2 



£1 

Pi 



B(ai; 1 - Ai)(V>(l + ai - Ai) - V(«i)) 



C2 

Pi 



B(a 2 ; 1 - A 2 )(^(l + a 2 - A 2 ) - ip(a 2 )) - p. 



If Ai or A 2 G {1,2} i/ie characteristic exponent can be computed using the 
following two integrals: 



(41) 



(e ixy - 1 - ixy) 



-a/3x 



1 - e~P x 



dx 



02^ («) 
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— afix 



(42) 



(e txy - 1 - ixy) 



o 



(1- 

— 1 - a + — 







£2 ^ («)■ 



Proof. First we assume that A G (0, 1). Performing change of variables 
ti = exp(— and using integral representation for beta function (formula 
8.380.1 in [20]) we find that 







1 — ixz) 



(1 



dx 



iz 



a--;l-A 



B(a;l - A) 



dz 



-^B( a- ^;1-A 



/3 



z=0 



and we obtain the desired result (40). The left-hand side of the above equa- 
tion is analytic in A for Re(A) < 3, and the right-hand side is analytic and 
well defined for Re(A) < 3, A 7^ {1,2}; thus by analytic continuation they 
should be equal for A £ (0, 3) \ {1, 2}. 

Assume that A = 2. Then using binomial series we can expand 

(1 — exp(— x))~ 2 = ^^{n + 1) exp(— nx), 

n>0 

which converges uniformly on (e,oo) and obtain 



/3 



1 — ixz)- 



■ dx 



E 



n>0 



(l_ e -/3x)2 

n + 1 n + 1 iy n + 1 



n + a — iy/(3 n + a (3 (re + a) 2 



n>0 



1 



1 



fa-«)E 



n + a — iy/ (3 n + a 
1 



/3 



n>0 



(n + a) 



and using the series representation for digamma function (formula 8.362.1 
in [20]) we obtain (42). Derivation of formula (41) corresponding to the case 
A = 1 is identical. □ 



The following theorem is the analogue of Theorem 5, and it is the main 
result in this section. 



Theorem 10. Assume that q>0 and that *S(z) is given by (40)- 
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(i) Equation q + ^f(iQ = has infinitely many solutions, all of which 
are real and simple. They are located as follows: 

C " E (-0101,0), 
C + e(0,/3 2 a 2 ), 

(43) 

Cn E (0 2 (a 2 + n- l),0 2 (a 2 +n)), n > 1, 

Cn E (01 (-ai + n),0i(-ai +n + 1)), n < -1. 

(ii) If o~ j^O we have 

Cn+i = 2 (n + a 2 ) + - (n + a 2 ) A2 ~ 3 

cr 2 0|r(A 2 ) 



(44) 



+ 0(n A2 "^ £ ), ra^+oo, 
2ci 

Cn-i = 01 (n - an) - ^2g2f()^) (~ n + ai ) Al_3 



+ 0(n Al_3_e ), n->--oo. 

(iii) If a = we have 

( n+ s = 02 (n + a 2 + w ) + ^4(n + a 2 + w ) A 

(45) 

+ 0(n A e ), n-^+oo, 

where coefficients wo and A are presented in Table 1, 5 £ {0, 1} depending 
on the signs of wo and A and 

i / . /ci^ra-Ai) 

xo = — arctan sin vrAo ; cosIttAt) 

The corresponding results for n — > — oo can be obtained by symmetry consid- 
erations. 

(iv) Function q(q + ^(z)) -1 can 5e factorized as follows: 

q 1 tt l + ^/(0 2 (n-l + a 2 )) 

q + V(z) (l + ^/C + )(l + ^/C ")4i 

(46) 

-q l + w/(0i(n + l-ai)) 



x 



where the infinite products converge uniformly on the compact subsets of the 
complex plane excluding zeros/poles of q + * (z) . 

Remark 5. When a = the remaining cases Ai < 2, A 2 = 2 and Ai = 2, 
< A 2 < 3 are not covered by Theorem 10. The interested reader can derive 
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these asymptotic expansions by using formulas (41), (42) and the following 
results for the digamma function (see formulas 6.3.7 and 6.3.18 in [1]): 

■0(1 - Z) = tjj(z) + 7TCOt(7T,z), lj)(z) = lll(z) - — + 0(z~ 2 ), 



2z 



oo. 



Proof of Theorem 10. The proof of (i) is very similar to the corre- 
sponding part of the proof of Theorem 5. We separate equation q + ^f(iC) = 
into a jump part and a diffusion part, find points where the jump part goes 
to infinity and by analyzing the signs we conclude that on every interval 
between these points there should exist a solution. 

The proof of (ii) and (iii) is based on the following two asymptotic for- 
mulas as C — > +oo: 

T(2o + 7-l) +0(c - 2) 



B(a + C;7)=r( 7 )C" 



2C 



B(«-C;7)=r (7) si ll c - a -^ c^ 



i + 7 (2a + 7 -l) +Q(r8) 



sin(7r(C — a)) 

The first asymptotic expansion follows from the definition of beta function 
(39) and formula 6.1.47 in [1], while the second formula can be reduced to 
the first one by applying a reflection formula for the gamma function. 

If a 7^ and £ — > +oo we use (40) and the above formulas and rewrite 
equation q + = as 



sin(7r(C//32 - «2 + A 2 )) 
sm(n((/l3 2 - a 2 )) 



C 3_A2 +o(c i_A2 ) 



while if a = and £ 

sin(7r(C//32 



a 2 



2c 2 r(l - A 2 ) 

+ o(c Al_A2 ) 

+oo we have 
A 2 )) _ 



1-A 2 n 



>2-A 2 + ci/3 2 2 r(l - Ai) Al _ Aa 
c 2 ^V{\ - A 2 ) 



sin(^(C//3 2 - a 2 )) c 2 r(l - A 2 ) 



Table 1 

Coefficients for asymptotic expansion of £ n when a = 



co 



Ai<2, A 2 <2 
Ai < 2, A 2 > 2 
Ai > 2, A 2 < Ai 
Ai > 2, A 2 > Ai 
Ai>2, A 2 = Ai 





2- A 2 


2- A 2 
x 



sin(-n-A 2 ),g|p 
" 7rc 2 r(l-A 2 ) 



cj/Sj 1 T(l-A!)r(A 2 ) 
sin(^A 2 ) ei/33 1 + 1 r(l-Ai) 
c 2 /3^r(l-A 2 ) 

-P^f^f F(A 2 ) 



A 2 -2 
2- A 2 
A 2 — Ai 
Ai — A 2 
2- A 2 
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+ o(c 1 ~ Xa ) + o(c Xl ~ Xa ~ 1 ). 

Asymptotic expansions (44) and (45) can be derived from the above formulas 
using the same method as in the proof of Theorem 5. 

In order to prove factorization identity (46) and the fact that there are no 
other roots, we use exactly the same approach as in the proof of Theorem 
5. Again we choose an entire function P(z) which has zeros at the poles of 
ty(z) with the same multiplicity, and the choice is obvious due to (40): 



(47) P{z) 



,,ai -f) r r + ! 



and function Q(z) is defined by as (? _1 (1 + ^(z))P(z). Function P(z) can 
be expanded in infinite product using Euler's formula (see formula 6.1.3 in 
[1]). Using (40) and (47) and asymptotics for gamma function we find that 
function Q(z) has order equal to one, and thus again we can use Hadamard's 
theorem to expand it as infinite product, and finally we use asymptotics for 
infinite products supplied by Lemma 6 and asymptotics for £ n given by (44) 
and (44) to prove factorization identity (46). □ 

We can also derive a result similar to Lemma 7 using the entire function 
Q(z) defined by (47). While there is no closed form expression for derivatives 
of gamma function, they can be easily computed numerically. Our final result 
in this section is the analogue of Theorem 8, and the proof is identical. 

Theorem 1 1 . For q > 

l + iz/(p2(n-l + a 2 )) 



^ W l + l + iz/C n 

+ _ 1 rr 1 +iz/(/3i(n + 1 -ai)) 

Infinite products converge uniformly on compact subsets o/C\zR. The den- 
sity of S r is given by 

d : nSr <x) = -%Co^ X - E c kCke^ x , 



dx 



fe<-i 



where 



n 



l-C -/(/?i(n + l-ai)) 



1-C "/Cn 



n<-l 



l-Cfc/C3i(A; + l-ai)) TT l-(k/((3i(n + l- ai )) 

1 - c^/Co- n Ji^ k i-a/Cn 
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Remark 6. In the case a = 0, Aj < 2, and p > we have a process of 
bounded variation and negative drift, and thus the distribution of S T will 
have an atom at zero, which can be computed using the following formula: 

P(5 T = 0)= Urn E[e~ zST }= lim </>+(iz) 

= TT 2H 

ai @ l n <-l ~ 

Using asymptotic relation (45) one can see that the above infinite product 
converges to a number between zero and one. 

5. Implementation and numerical results. In this section we discuss im- 
plementation details for computing the probability density function of S T r q \ 
and St- In order to illustrate the main ideas we will use the process X% 
defined in Section 3; however, the implementation for a general Xt from 
the /3-family would be quite similar. Our main tools are Theorem 8 and 
asymptotic expansion for ( n given in Theorem 5. 

First let us discuss the computation of density of S T . The first step would 
be to compute solutions to equation q + ^(iQ = 0, and for q real this is 
a simple task: for n large we use Newton's method which is started from 
the point given by asymptotic expansion (16) or (17). To compute or ( n 
with n small we use localization result (15) and the secant (or bisection) 
method to get the starting point for Newton's iteration. Overall this part of 
the algorithm is very computationally efficient and can be made even faster 
if we compute different £ n in parallel. 

The second step is to compute coefficient c^7 which are given by (37). Each 
term in the infinite product is 1 + 0(n~ 2 ); however, as we show in Propo- 
sition 12, we can considerably improve convergence by using our knowledge 
of the asymptotic expansion for Q n . The final step is to compute the density 
of S T using formula (36). Note that the series converges exponentially for 
x > 0. When x is small the convergence is slow, and the asymptotic behavior 
as x — > + would depend on the decay rate of coefficient ; however, we 
were unable to prove any results in this direction. 

Proposition 12. Assume that Q n = n + (3+ + u+p)* + 0(n~ 3 ) as 
n — > +00. Then as N — > +00 we have 



n 

n>N 



(48) 



1 + zj(n + a) 

1 + Z/Cn 

T(N + a)T(N + f3 + z) 
T(N + P)T(N + a + z) 
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x exp[A! (/^(ft ft N) - f 1A (z + ft ft iV)) 

+ A 2 (f 1 , 2 (p, ft iV) - / 1)2 (z + ft ft iV)) + 0(N~ 3 )], 
where f ai ,a 2 ( z ij z 2', N) can be computed as follows : 

fa 1 ,a 2 (zi,Z 2 ;N) 

ai + a 2 + k-l ' 

(49) 

1 ai ao 

77 + T777 777 + 



2 12(zi+JV) 12(z 2 + iV) 



+ (zi + Af)^ ai (^ + A r )^ 2 
+ 0(N- ai - a2 ~ 3 ). 
Proof. First we define for ct\ + a 2 > 1 

n>7V 

The proof of the asymptotic expansion (49) is based on the Euler-Maclaurin 
formula 

n>N JN 

where we take /(x) = (x + zi)~ Ql (x + 22)"°+ To obtain (49) we compute 
the integral by changing variables y = {x + zi) -1 , expanding the resulting 
integrand in Taylor's series at y = and integrating term by term. 

In order to obtain formula (48) we follow the steps of the proof of Lemma 

6 

-pr 1 + z/(n + a) _ -i-r l + z/(n + a) -i-r l + z/(n + 0) 

_ r(AT + a)r(A^ + /3 + z) -p. l + + 

"r(iv + /3)r(iv + a + z) n ll i + z/Cn ' 

Next, we approximate ln(l + uj) = lo + 0(w 2 ) as u — > and obtain 

I>(W) 

nt^ V" + /V 4^ + * + 
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n>N 



A x A 2 _ 4 , 

+ 7Z i ay, + °( n . 



(n + /3) 2 (n + P) 3 



E 



+ 7 ^7 ^ + 0(n" 



(z + n + /3)(n + /3) (z + n + /3)(n + (3f 



n>N 

which completes the proof. □ 

Computing the density pt(x) = 33: IP (St < x) of St at a deterministic time 
t > requires more work. Our starting point is the the fact that the density 
of S r [which we denote by p (q,x)] is the Laplace transform of q x pt{x) 

d d r°° 

p S (q, x) = ^nS T{q) < x ) = —j^ qe -*V(St < x) dt 



= q / e- qt p t {x)dt. 
Jo 

Thus pt(x) can be recovered as the following cosine transform: 

p S (qo +iu,x) 



2 Z -00 
Pt ( x ) = -e<?°* / Re 

7T Jo 



qo + iu 



cos(tn) du, qo > 0. 



We see that to compute this Fourier integral numerically we need to be able 
to compute p (q, x) for q lying in some interval gG [qo, qo + iuo] in the com- 
plex plane. The main problem is that we need to solve many equations (22) 
with complex q. While the asymptotic expansions for Q n presented in The- 
orem 5 are still true, we do not have any localization results in the complex 
plane. It is certainly possible to compute the roots using the argument prin- 
ciple, and originally all the computations were done by the author using this 
method. However as we will see, there is a much more efficient algorithm. 

We need to compute the solutions of equation q + ^(iC) = f° r an Q G 
[qo,Qo + iuo\. First we compute the initial values: the roots , ( n for real 
value of q = qo using the method discussed above. Next we consider each 
root as an implicit function of u : Cn{ u ) is defined as 

*«„(«)) + (qo + m) = 0, Cn(0) = Cn- 
Using implicit differentiation we obtain a first order differential equation 

d(n(u) _ 1 

du ~ tf'(iC„(u)) 

with initial condition Cn(0) = Cn- We compute the solution to this ODE using 
a numerical scheme, for example, an adaptive Runge-Kutta method, and at 
each step we correct the solution by applying several iterations of Newton's 
method. Again, for different n we can compute Cn(u) in parallel. 
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Figure 2 shows the result of this procedure. We have used the following 
values of parameters: a = 1, fi = —0.1 and a = 0.25 and computed zeros 
and Cn for q € [1,1 + 200i]. The graph shows interesting qualitative behavior: 
all zeros except Co converge to the closest pole of at a + n, while 

has no pole nearby [since is regular at £ = a] and it converges to oo 

while always staying in C + . If a < the situation is exactly the same, except 
that now escapes to oo while always staying in C~. We have repeated 
this procedure for many different values of parameters, and from this numer- 
ical evidence we can make some observations/conjectures. It appears to be 
true that the roots never collide which means that we have no higher order 
solutions to q + ^(z) = for all q € C. It also seems that the roots never 
cross the real line. All these observations are based on numerical evidence, 
and we did not pursue this any further to obtain rigorous proofs. However 
there is one fact that we can prove rigorously: there are not going to appear 
any new, unaccounted zeros. This could be proved by an argument that we 
have used in the proof of Theorem 5 to show that there are no extra zeros. 

The results of our computations are presented in Figure 3. The parameters 
are a = 1, fj, = —0.1 and a = 0.25; the surface p s (q, x) is on the left and pt(x) 
on the right. 

6. Conclusion. In this paper we have introduced a ten-parameter family 
of Levy processes characterized by the fact that the characteristic exponent 
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is a meromorphic function expressed in terms of beta and digamma func- 
tions. This family is quite rich, and, in particular, it includes processes with 
the complete range of behavior of small positive/negative jumps. We have 
presented results of the Wiener-Hopf factorization for these processes, in- 
cluding semi-explicit formulas for Wiener-Hopf factors and the density of 
the supremum process S T . 

These Levy processes might be used for modeling purposes whenever one 
needs to compute distributions related to such functionals as the first passage 
time, overshoot, extrema, last time before achieving extrema, etc. Some 
possible applications in Mathematical Finance and Insurance Mathematics 
include pricing barrier, lookback and perpetual American options, building 
structural models with jumps in Credit Risk, computing ruin probabilities, 
etc. 

Finally, we would like to mention that one can use the methods presented 
in this paper, and, in particular, the technique to solve q + ^/(z) = for 
complex values of q discussed in Section 5, to compute Wiener-Hopf factors 
arising from a more general factorization identity (see Theorem 6.16 in [21]) 

$+(w,z)t?~(w,z), w,zeR, 



q-iw + ^(z) ~ qK ' 1 q 
where 

¥+(«/, z) =E[e iw ^ +izS % V-(w,z) = E[e iw ^ +izI T], 

and Gt (G t ) are defined as the last time before t when maximum (minimum) 
was achieved 

Gt = sup{0 <u<t: X u = S u }, G t = sup{0 <u<t: X u = I u }. 
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